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ABSTRACT 

On 2011 May 30, quasi-periodic fast propagating (QFP) magnetosonic waves 
accompanied by a C2.8 flare were directly imaged by the Atomospheric Imaging 
Assembly instrument on board the Solar Dynamics Observatory. The QFP waves 
successively emanated from the flare kernel, they propagated along a cluster of 
open coronal loops with a phase speed of ^834 km s~^ during the flare's rising 
phase, and the multiple arc-shaped wave trains can be fltted with a series of 
concentric circles. We generate the k-uj diagram of the Fourier power and flnd 
a straight ridge that represents the dispersion relation of the waves. Along the 
ridge, we flnd a lot of prominent nodes which represent the available frequencies 
of the QFP waves. On the other hand, the frequencies of the flare are also 
obtained by analyzing the flare light curves using the wavelet technique. The 
results indicate that almost all the main frequencies of the flare are consistent 
with those of the QFP waves. This suggests that the flare and the QFP waves 
were possibly excited by a common physical origin. On the other hand, a few low 
frequencies revealed by the k-u diagram can not be found in the accompanying 
flare. We propose that these low frequencies were possibly due to the leakage of 
the pressure-driven p-mode oscillations from the photosphere into the low corona, 
which should be a noticeable mechanism for driving the QFP waves observed in 
the corona. 
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INTRODUCTION 



The magnetically dominated plasma of the solar corona, which is an elastic and com- 
pressible medium, can support the propagation of various magnetohydrodynamic (MHD) 
waves (e.g., slow mode, Alfven, and fast mode waves). In the last decade, observations from 
ground-based and sp ace-borne instruments have led to the detection of various waves in the 



solar atmosphere (see | De Moortel et al.ll2002bl : iNakariakov fc Verwichtdl2005l : iBanerjee et al. 



20071 : IPe Moortell l2009l ) . For example, (1) oscillat ions or standing waves o 



modes in a variety of coronal magnetic structures (lAschwanden et al 



19991 : IWang et al.ll2002l : INakariakov et al.ll2003l : lOkamoto et al.ll2007l : lOfman fc Wang 



1999 



slow and fast 



Schriiver et al. 



(2) propagating compressible s low magnetosonic wa ves in solar polar plunies (IDeForest fc Gurman 



1998 : 



Ofman et al.lll997 



20001 . l2002al : 



Marsh et al 



2003 



Robbrecht et all 



19991) and in coronal loops (IBerghmans fc Clettdl 19991 : 



incompressible Alfven waves in the corona (ICirtain 



20011: 



Wang 



20071 : iTomczvk fc Mcintosh 



waves in coronal lo ops ( 1 Williams et al. 



tions in 



2009; 



Jess et al 



2001 



et al.ll200 



et al. 



200 



2009a. 



Km 



De Moortel et al 



transverse 



De Pontieu et al.ll2007l : ITomczvk et al 



20091) , and (4) compressible fast magnetosonic 



20021 : McLaughlin et al.ll2009l : iLiu et all 12011 



Ofman et al.ll2011r>. It should be noted that the Alfven wave interpretation of the observa- 



Tomczyk et al.l ( 120071 ) could al ternatively be explaine d as kink waves ( IVan Doorsselaere et al. 



20081). The single-pulse "EIT waves" (IThompson et al.lll998h with tvpical speeds of 200-400 



km s ^ were often interpr e ted as fast mag netosonic waves (IWu et al 



200ll : lOfman fc Wang 



2OO2I : IVeronig et al.l I2OIOI : iLiu et al.l l2010r) and th ought to be the coronal counterparts of 



the chromospheric Moreton waves jMoretonl[l960h . which have be en successfully explained 



as coronal shock waves sweeping the chromosphere ( lUchidal Il968l ). However, the wave in- 



some drawbacks in 


the wave models (see 


Delannee & Aulanier 


1999 




Delannee et al. 


2008; 


Klassen et al. 


2000 




AttriU et al. 


2007, 


2009: 


Wills-Davevlbood: 


WarmuthI 


2OI0I; 


Chen et al. 



2002 



200,4 



Chen fc Wu 



201 ll ). Very recently, the universal wave interpretations for coro- 
nal oscillations are challenged by the discovery of the rapid quas i -periodic mass up-flows 



at transition- region and coronal temperatures (|Sakao et al 



Tian et al 



De Pontieu fc Mcintosh 2010; Mcintosh fc De Pontieu 2009a 



2007 



De Pontieu et al. 



Mcintosh et al 



2010. 



2009; 



2011 



_ ^ . I 17 , — = , I ■>^7 I 1 I 17 I 17 

201 ll ). and those up-flows were suggested to provide hot plasm a into the corona 



and thereby heating the corona or as a poss i ble source for the fast solar wind f|De Pontieu et al 



20091. 2011 



201( 



Mcliitosh fc De Pontieull2009bl : lPetedl2010l : iDe Pontieu fc Mcintosh 
Tian et al.ll2"oill ). 



201 



3 



Hansteen et al 



On the other hand, the oscillations and quasi-periodic pulsations (QPPs) of solar flares 
have also been detected for many years in wide wavelength bands rang ing from radio to 
gamma-ray, on timescales from fractions of seconds to several minutes ( e.g.. iParks fc Winckler 



19691 : lAschwanden I Il987l : iFleishman et al.l l2002l : iFouUon et al.l l2005l : iKislvakov et al. 



2006 
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Tanl 120081: i Nakariakov fc Melnikovl I2OO9I : iNakariakov et al.l I2OI0I : iKuprivanova et all I2OI0I : 
Tan fc Tanl[20l3 ). QPPs can usually be observed in light curves of solar flares in various 
wavelength bands during the entire flaring pro cess. Moreover, they can also b e observed in 
stellar flares with a period of several minutes ( IMathioudakis et al.l 120031 . l2006l ). In different 



wave l ength bands , QPP s often shows a synchronous evolvement patten in phase (lAsai et al. 



200 ll : llnglis et al.l l2008l ) . In som e events, several signiflcant periods could be detected si- 
multaneously i n one single flare (iKislvakov et al.ll2006 l). For different timescales of QPPs 
in solar flares, Nakariakov fc Melnikov ( 20091 ) split them into short (sub-second) and long 
period QPPs. The short period QPPs are usually observed in radio emission, and are likely 
associated with the interaction of elect ro-magnetic, plasma or whistler waves with acceler- 
ated particles (e.g., lAschwanden Ill987l ). while the long period QPPs are usually detected in 
microwave , white light, EUV and X-ray e missions, and are p ossib ly associated with MHD 
processes (INakariakov fc Melnikovl |2009|) . INakariakov et al.l ( 120041 ) performed a numerical 
study on the evolution of a coronal loop in response to an impulsive energy release, flnding 
that the QPPs of the loop density is associated with the second ha rmonic of a standing 
slow magnetosonic wave. In another publication ( ITsiklauri et al.ll2004l ). it is further pointed 
out that the generation of the QPPs of the loop density is independent of the location 
of the heat deposition in the loop. Although there are many observational and theoreti- 
cal studies in the past several decades, the underlying physical mechanisms that generate 
QPPs in solar flares still remain an open question. Up to the present, there are several 
possible inte rpretations for this issue, including modulation of electron dynamics by MHD 



oscillations ( jZaitsev fc Stepanov 



relea s es by MHD waves (iRoberts et al 



1982 



20061 : llnglis et al.ll2008l : lSvch et al.ll2009l ). MH D flow overstabilitie s ( 



1198 



1984: 



Fleishman et al.ll2008h. periodic tr i ggering of energy 



FouUon et al. 



oscillatory regimes of magnetic reconnections (I Kliem et al.ll2000l : iMcLaughlin et al.ll2012l ). 



2005: 



Nakariakov et al. 



Qfman fc Sui 2006.). a nd 



Quasi-periodic propagating fast magnetosonic waves with high speed remain the least 
observed among all coronal MHD waves. The scarcity of fast wave observations was mainly 
due to the instrumental limitation in the past. For wavelengths that are comparable with 
the characteristic size of coronal magnetic structures (e.g., length of coronal loops, width 
of coronal polar plumes), the typical periods of the wa ves are in the range from a few 
seconds to several minutes ( Nakariakov fc Verwichte 20051 ) . Hence both spatial and temporal 
resolutions are necessary ingredients for the detection of fast magnetosonic waves. Namely, 
the time cadence must be shorter than the periods and the pixel size must be smaller than 
the wavelength of the fast magnetosonic waves. The flrst imaging observation that had 
been interpreted in terms of propagating fast magnetosonic wave had been reported by 
Williams et al.l ( 120011 . 120021 ). They successfully imaged a propagating fast magnetosonic 
wave during the total solar eclipse on 1999 August 11, in which the wave travels through the 
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apex of an active r egion coronal loop with a speed of 2100 km s ^. ICooper et al.l ( 120031 ) and 
Nakariakov et al.l ( 120031 ) modeled the evolution of the propagating fast mag netosonic wave 



and c o nfirni ed the formatio n of the quasi-pe riodic wave trains predicted by [Roberts et al. 



(119831 . Il984j ). Very recently, iLiu et al.l ( 120111 ) reported the quasi-periodic fast propagating 
(QFP) magnetosonic waves in the low corona usi nR the EUV obser vations taken by the 
Atmospheric Imaging Assembly instrument (AIA; iLemen et al.l l2012l ) on board the Solar 
Dynamics Observatory (SDO). The authors found that the multiple arc-shaped wave trains 
emanated consecutively near the associated flare kernel and propagated outward along a 
funnel of coronal loops with a phase speed of ^2200 km s~^. In addition, they also found 
that the main period of the associated QPPs flare was almost the same with the strongest 
period (181 s (5.5 mHz)) of the QFP waves. They thus concluded that the flare and the 
QFP waves have a common origin of possibly quasi-periodic magnetic r econnectio n . A th ree- 
dimensional numer ical simulation stud y of the QFP waves reported by iLiu et al.l ( 120111 ) has 
been presented by lOfman et al.l (120111 ) , in which they successfully reproduced the multiple 
arc-shaped wav e trains that ha ve similar amplitude, wavelength, and propagation speeds 
as reported bv iLiu et al.l ( 12011 ). They also discussed the possible excitation mechanism 
of the QFP waves and the applications of the obser vations for cor o nal s ei smology. Other 



(2009 


), 


Fedun et al. 


(2011) 



In this paper, we present an observational study of a QFP wave event that occurred on 
2011 May 30. It was accompanied by a GOES C2.8 flare in NOAA active region AR11227 
(S19E57) and a faint coronal mass ejection (CME) observed in white light. We flnd that 
the QFP wave trains consecutively emanated from the flare kernel and propagated along a 
bunch of open active region coronal loops which extended from AR11227 to the heliosphere, 
while another propagating wave are observed in the nearby closed transequatorial loops 
that connect the neighbouring active region AR11226 and the outskirts of AR11228 in the 
northern hemisphere. Using high spatial and temporal observations taken by AIA, we flrst 
generate the k-uj diagram of the QFP waves with Fourier transform, and then isolate the 
periods of the associated pulsation flare with light curves in various wavelength bands. By 
comparing the periods (frequencies) of the wave and the flare, we try to obtain some new clues 
about the driving mechanism of the QFP magnetosonic waves. Observations are described in 
Section 2, results are presented in Section 3, discussions and conclusions are given in Section 
4. 
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2. OBSERVATIONS 

The AIA on board SDO is propitious to detect fast propagating features such as fast 
magnetosonic waves of short periods. The telescope captures image of the Sun's atmosphere 
out to 1.3 i?0 and has high time resolution of up to 12 s. AIA has high spatial resolution 
of 1" .2 in seven EUV and three UV-visible channels, which cover a wide temperature range 
from logT = 3.7 to logT = 7.3. All these parameters of AIA are necessary ingredients for 
detecting fast propagating waves in the low corona. In the event analyzed here, the QFP 
waves along open active region coronal loops can only be observed in the AIA 171 A (Fe IX; 
logT = 5.8) observations, which mainly image the quiet corona and the upper transition 
region. In the meantime, another wave that was trapped in closed transequatorial loops 
can be detected at both 171 A and 193 A (Fe XII : IorT = 6.1) l ines. The Reuven Ramaty 



High Energy Solar Spectroscopic Imager ( RHESSL iLm et al.ll2002l ) hard X-ray (HXR) count 
rates in the energy bands (4 s integration) of 12-25 keV and the Geostationary Operational 
Environmental Satellite (GOES) soft X-ray (SXR) flux are also used to analyze the periods 
of the accompanying QPPs flare. All images used in this paper are differentially rotated to 
a reference time (11:00:00 UT), and solar north is up, west to the right. 



3. RESULTS 

The event presented in this paper includes the QFP waves along the open coronal loops 
that rooted in AR11227, and another fast magnetosonic waves in the closed transequatorial 
loops that connect the neighbouring active region AR11226 (S21E42) and the outskirts of 
AR11228 (N17E55) in the northern hemisphere. These waves were accompanied by a GOES 
C2.8 flare that occurred in AR11227 on 2011 May 30. The start, peak, and end times of the 
flare are 10:48, 10:57, and 10:59 UT respectively. According to the RHESSI &.diTe record, the 
heliographic coordinates of the flare kernel were x = —787'^ y = — 283'^ In addition, a faint 
CME in white light was also observed to be associated with a jet-like plasma eruption which 
occurred in AR11227 just before the flare and the waves. 

An overview of the event is displayed in Figure [H One can see that a cluster of funnel- 
like open loops rooted in AR11227 and extended northeast off the disk limb, and another 
bunch of closed transequatorial loops connected AR11226 and the outskirts of AR11228 in 
the northern hemisphere. These loops can be identifled clearly on the raw and the flltered 
AIA 171 A images (see Figure [l](a) and (b)). The QFP waves were observed along the open 
loops. They can be distinguished at about 10:50:00 UT and then multiple arc-shaped wave 
trains are noticed from the 171 A running difference images (see Figure [l](c) and (d), and 
Animation 1 available in the online version of the journal). The wave trains consecutively and 
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alternately emanated from the flare kernel and exhibited an arc-hke shape. They propagated 
northeastward along the open loops and flnally dimmed off the disk limb successively. This 
suggests that the QFP waves were not propagating; surf ace waves as the Moreton and EIT 



pre 



waves, in agreement with the results bv iLiu et al.l ( 120111 ). The start time of the QFP waves 



(10:50 UT) was slightly lagging the beginning of the accompanying flare (10:48 UT), and 
the end time of the QFP waves (11:02 UT) was also lagging the flare peak time (10:57 UT) 
several minutes. This close temporal relationship between the QFP waves and the flare 
implies that they might be closely connected. In addition, we note that the QFP waves 
could only be observed in AIA 171 A observations, and the open loops that guided the QFP 
waves could not be observed in the other wavelength bands. These results indicate the subtle 
temperature dependence of the QFP waves. 

It should be noted that a jet-like plasma ejection around the base of the open loops was 
observed before the beginning of the start of the accompanying flare and the QFP waves. This 
small plasma ejection started at about 10:30 UT and ended at about 10:53 UT. It erupted 
outward eastwardly and it possibly resulted in the faint CME due to the intimate temporal 
and spatial relationship between them (see the long white arrow in Figured]). Furthermore, 
the obvious coronal dimming region close to the source region of the plasma ejection was 
possibly the on-disk signature of the occurrence of the CME. The discrepancies on the start 
times, the propagating directions of the QFP waves and the jet-like plasma eruption together 
indicate that the two phenomena did not have a cause and eflFect relationship. In this paper, 
we mainly focus on the relationship between the QFP waves and the associated flare. 

For the multiple arc-shaped wave trains propagating along the open loops, each of them 
can be fltted with a circle whose center is located at the site of the flare kernel. To analyze the 
kinematics of the QFP waves, we place flve lines (A0-A4) in a sector region where the QFP 
waves were propagating. The angle between the adjacent lines is 10° (see Figure [T](c)). By 
averaging 10 pixels across each line and then stacking the obtained proflles in time sequence, 
the time-distance diagrams along the line can thus be obtained, and the results are shown 
in Figure [21 For the wave guided by the closed loops, it propagated resembling a surface 
wave that was conflned on the solar surface(see Animation 2). In addition, obvious coronal 
dimming region can be identifled following the wave fronts (see Figure [l](e) and (f)). To 
minimize human subjectivity and the spherical projection effect, we select the flare kernel 
as the new "north pole" and the angle between the new heliographic "longitudes" is still 10° 
(see Figure [11(e)). Along each "longitudinal" great circle, we obtain the intensity proflle of 
each image by averaging 10 pixels in the "latitudinal" direction. By stacking these intensity 
proflles over time, one can obtain a two-dimensional time-distance diagram, as shown in 
Figure [71 
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Figure [2] shows the time-distance diagrams obtained from AIA 171 A running (Fig- 
ure [2](a)-(c)) and base (Figure [2](d)-(f)) difference images along cuts A1-A3 shown in Fig- 
ure [11(c). These time-distance diagrams exhibit the kinematics of the QFP waves guided by 
the open coronal loops. From the running difference time-distance diagrams, one can see 
that the multiple wave trains are shown as a series of steep stripes with positive slope repre- 
senting the propagating speed of the QFP waves. It can be seen that these steep stripes are 
originated from the flare kernel and become prominent at a distance of about 100 Mm far 
from the flare kernel. It should be noted that the QFP waves underwent a deceleration phase 
before t hey became promin ent. Since the solar corona is a highly inhomogeneous medium in 



nature (jOfman et al.ll201ll ). the deceleration of the QFP waves was possibly resulted from 
the changing of the properties of the propagation medium (such as the variation of density 
and magnetic field strength along with height of the loops) . The speed of each stripe in the 
plane of the sky is measured using a linear fit method. Within the distance ranging from 120 
to 220 Mm far from the flare kernel, within which the wave trains were mostly prominent, 
we find that the speeds of the different stripes show little difference along different lines, and 
the average speed of all wave trains is about ^774 km s~^. For the wave trains within the 
distance less than 100 Mm far from the flare kernel, the measurement is inaccurate since the 
wave trains are very obscure. For a rough estimation, the speed of the wave trains within 
such a distance should be higher than 1000 km s~^. It should be noted that a dimming 
region is detected on the base difference time-distance diagrams close to the flare kernel (see 
the bottom row in Figure [J]). As mentioned above, it was associated with the jet-like plasma 
ejection and the CME observed in white light, and it had probably no direct relation with 
the QFP waves studied here. 

We further check the intensity variation of the QFP wave trains along and cross the lines 
that are used to obtain the time-distance diagrams, and the results are shown in Figure [3l 
We cut a narrow subregion from the running difference time-distance diagram along A2 
(Figure [3](a)). The time interval for this subregion is from 10:56:36 UT to 10:58:00 UT, 
corresponding to seven time points of the AIA's 12 s time cadence. On this sub time- 
distance diagram, the wave trains are more clear. By making a linear fit to the wave trains, 
we obtain the speed of the propagating wave which is ^724 km s~^. The intensity profile at 
each time point as a function of distance shows obvious quasi-periodic variation as indicated 
by the red curves in Figure [3](b) , in which each profile is normalized into [-6,6] DN and 
then incrementally shifting by 12 DN that corresponds to the time cadence of the AIA 
instruments. Thus the x-axis also servers as the elapsed time. It is interesting that each 
intensity profile can be fitted with a sine function Asin[27T{r — ro)/X] (blue curve), here A 
being the amplitude, A the wavelength, and Tq the initial phase in distance. The average 
fitted parameters are A = 23.8 Mm, A = 2.58 DN. Based on these fitted curves, we can 
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obtain the phase speed of the QFP waves in the plane of the sky, which is ^723 km s"-'^, 
almost the same with the phase speed measured from the sub time-distance diagram shown 
in Figure [3](a). To investigate the intensity variation across the lines, a subregion is cut 
from the base difference 171 A time-distance diagram obtained from A2 (dashed box shown 
in Figure [2](e)), and the intensity profiles at each position is plotted as a function of time, 
in which the next profile of larger distance is stacked on the previous one. The QFP wave 
trains are more clear, and obvious periodic variation of intensity can be found (see the black 
arrows in Figure [3](d)). 

For MHD waves generated impulsively, we often observe the propagating wave packets 
in imaging observations. These wave packets represent a Fourier integral over all frequencies. 
Therefore, one can decompose these frequencies with the method of Fourier transform. In 
line with this thought, we first extract a three-dimensional data cube in the time domain 
(i.e., in (x,^,t) coordinates, the time interval is from 10:48^11:00 UT, corresponding to 
the rising phase of the associated fiare), then transform the data cube into the frequency 
domain (i.e., in (fcx, fcy, ly) coordinates, where and ky are wavenumbers along x and ^, u is 
frequency). Then, the Fourier power is summed in the azimuthal 9 direction of cylindrical 



DeForest 


2004; 


Liu et al. 


2011) 



-id 



diagram is Ak = 6.83 x 10~^ Mm~^, Au = 1.39 mHz, which are calculated based on the FOV 
of the selected region and the time interval of the data cube. An obvious liner steep ridge 
that represents the dispersion relation of the QFP waves can be found in the k^u diagram 
(see Figure [4](a)). This ridge can be well fitted with a straight line passing through the 
origin up to the Nyquist frequency of 41.7 mHz. This gives the phase ('Uph = ly/k) and the 
group {vgj. = dv/dk) velocities of the QFP waves (^834 km s~-^). It should be noted that 
many prominent nodes can be found on this bright ridge, and they represent the available 
frequencies (periods) of the QFP waves. We also plot the normalized intensity profile of 
this bright ridge in Figure [4](a) as a red curve, from which we find the intensity peaks that 
represent the available frequencies (periods) of the QFP waves (see the horizontal dotted lines 
in Figure [4](a)). We find that the frequencies (periods) of the QFP waves range from 0.7- 
39.8 mHz (1428-25 s). Figure [4](b)-(d) are the snapshots of the Fourier filtered AIA 171 A 
running diflFerence images with a narrow Gaussian function centered at the peak frequencies 
of 11.2, 15.6, and 23.7 mHz, which highlight the corresponding QFP wave trains at different 
frequencies (also see Animation 3 available in the online version of the journal). Obviously, 
larger frequency corresponds smaller distance between adjacent wave trains, namely shorter 
wavelength A, since the propagating speeds of them are the same. 



Previous studies have indicated that an obvious source of propagating fast magnetosonic 
wave can be a fiare (IRoberts et al.lll984 ). So we measure the light curves of the associated 
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flare from the AIA observations by summing the intensities within a small region centered 
on the flare kernel. We show the sample light curves of 1600 A and 171 A in Figure O as well 
as the RHESSI HXR count rates in the energy band of 12-25 keV. The GOES 1-8 A SXR 
flux is also plotted in Figure [5](a). We flnd that all the light curves of the flare show obvious 
fluctuation and exhibit certain periodicity during the rising phase (10:48-11:00 UT) of the 
GOES C2.8 flare (Figure [5](b)). The periodic variation of the flare intensity become clearer 
by seeing the detrended intensity fluxes (Figure [5](c)), which are obtained by subtracting the 
corresponding smoothed fluxes using a 120 s boxcar. 

To isolate the periods of the accompanying pulsation flare from the detrended fluxes, we 
resort to the wavelet analysis method, a common technique for analyzing localized variations 
of power within a time series. This technique allows us to investigate the time dependence pe- 
riods within the observed data. The det ails of the procedure and the corresponding guidance 



are given by iTorrence fc Compd ( 119981 ). In our analysis, we choose the function "Morlet" 
as the basis function, and a red-noise signiflcance test is performed. Since both the time 
series and the wavelet function are flnite, the wavelet can be altered by edge effects at the 
end of the time series. The signiflcance of this edge effect is shown by a cone of influence 
(COI), deflned as the region where the wavelet power drops by a factor of e~^. Areas of the 
wavelet power spectrum outside the region bounded by the COI should not be included in 
the analysis. Our wavelet analysis results of the flare light curves are shown in Figure [6l in 
which we highlight the isolated frequencies (periods) of the QPP flare with horizontal red 
dashed lines. Note that the frequencies (periods) of the flare are determined by the global 
power where the signiflcance is above 95%. 

Based on the wavelet analysis technique, 16 different frequencies (periods) are isolated 
from the AIA and RHESSI light curves of the accompanying flare, and all of them are 
listed in Table [H It is interesting that a few periods (frequencies) are found simultaneously 
in different wavelength bands (e.g., 19 mHz (53 s)), while almost each wavelength band 
manifests multiple periods (frequencies) (e.g., 193 A). The period error of each period is 
determined by the full width at half maximum of each peak on the global power curve, 
which is obtained by fltting each peak with a Gaussian function. On the wavelet power 
spectrum diagrams, the blue contours represent the regions where the signiflcance is above 
95%. We measure the duration and cycles of each period deflned by these contours. The 
results show that the durations of these periods of the flare range from 166 s to 327 s, while 
the cycles range from 2.3 to 15.9. Since many measured periods are the same or showing 
little difference, we get the main periods by averaging the same or neighboring periods. This 
yields seven main periods (frequencies) of the flare (see Pfa (Ffa) in Tabled]). 

By comparing the periods (frequencies) of the pulsation flare and the periods (frequen- 
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cies) of the QFP waves, we find that five of the seven main fiare periods (frequencies) coincide 
weU with the QFP waves' periods (frequencies) which are revealed by the k-uj diagram shown 
in Figure [11(a). These corresponding periods (frequencies) are 25 ± 2.5 s (40 mHz), 32.6 ± 3 
s (31 mHz), 41 ± 3.8 s (24.4 mHz), 55.4 ± 5.8 s (18 mHz), and 83 ± 13.5 s (12 mHz). This 
result indicates that the pulsation fiare and the QFP waves are possibly excited by the same 
physical mechanism. On the other hand, the k^uj diagram of the QFP waves also reveals 
a few low frequencies (e.g., 9.3 mHz (108 s) and 2.5 mHz (400 s)) that can not be found 
from the pulsation fiare. We conjecture that these unmatched low frequencies are possibly 
caused by the leakage of the photospheric p-mode oscillations through the chromosphere and 
transition region into the low corona (e.g., 5 minutes oscillation), w hich has been considered 
as an obvious source of such low frequencies in a lot of studies (e.R., 



2002al : iMarsh et al.ll2003l : IPe Pontieull2004al . [20051; [Wang et al.[[2009a[ ]bl). 



De 



Moortel et al. 



2000. 



Apart from the QFP waves observed in the open loops, propagating waves along the 
closed transequatorial loops are also observed at both 171 A and 193 A wavelength bands. To 
minimize the infiuence from the neighbouring QFP waves that can only be observed on the 
AIA 171 A images, we investigate the propagating waves in the closed loops using only the 
AIA 193 A observations. The time-distance diagrams obtained from the AIA 193 A running 
difference images along the spherical curves that are shown in Figure [l](e) are displayed 
in Figure [3 From these time-distance diagrams, multiple steep stripes (representing the 
propagating wave trains) of both positive and negative slopes can be clearly observed on 
the time-distance diagrams (see Figure [71(b)). The speeds of all stripes in the plane of 
the sky along each cut are obtained by making a linear fit to the stripes. The average 
speed of all the stripes with positive slopes is ^240 km s~^, while it is ^220 km s~^ for 
the stripes with negative slopes which represent the refiection of the propagating waves at 
the northern footpoints of the closed transequatorial loops. In addition, a few wave trains 
changed their direction at a distance of about 200 Mm from the fiare kernel and their average 
speed is ^68 km s~^ (see the thin oblique arrow in Figure [71(c)). The refiection of the waves 
suggests that they should be real MHD waves rather than the so-called ^^pseudo- wave s" (e.g.. 



Delannee fc Aulanier 



19991 : 



Chen et al 



2002: 



AttriU et al.ll2007l : lDelannee et al.ll2008l ). which 



are suggested as the footprints or low coronal extensions of the associated CMEs. 



4. CONCLUSIONS AND DISCUSSIONS 

We investigate a QFP wave event in a great detail with the high temporal and spatial 
resolution observations taken by SDO/AIA^ in which the QFP waves showed multiple arc- 
shaped wave trains that were propagating along a cluster of open active region coronal loops 
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during the rising phase of the associated flare. We demonstrate that almost aU the periods 
(frequencies) of the accompanying pulsation flare coincide well with the periods (frequencies) 
revealed by the k-uj diagram of the QFP waves. This result highly suggests that the QFP 
waves and the accompanying pulsation flare were possibl y excited b y a c ommon physical 
origin. Similar result has also recently been proposed by Liu et al. ( 20111 ) in their report 
on the flrst QFP waves directly imaged by SDO/AIA. However, in their case, only one 
frequency of the QFP waves coincides with the accompanying flare. In our event, the perfect 
consistency of multiple frequencies of the pulsation flare and the QFP waves make us believe 
the existence of such a relationship between them, with a high confldence level. Furthermore, 
the k-u) diagram of the QFP waves also reveals a few low frequencies that can not be found 
in the flare, such as 9.3 mHz (108 s) and 2.5 mHz (400 s). We conjecture that these low 
frequencies are possibly the manifestations of the photospheric pressure-driven (p-mode) 
oscillations in the low corona. The leakage of the photospheric oscillations through the 
chromosphere and transition region into the low corona has been p r oved to be possib 



number of theoretical and observational studies ( Bel fc Lerov 



2002al: 



Marsh et al 



Didkovskv et al.ll201lh . 



2003; De Pontieu 2004a. 2005 



1977: 



Chen fc Priest 



De Moortel et al 



20061 : IWang et al.ll2009al ]b 



e m a 



2000 



It should be noted that the projection effect may affect our analysis results about the 
waves in the closed transequatorial loops to a certain extent. Since these loops were located 
within the solar limb, we can not distinguish the height of the propagating waves based on 
a single view point observations, since the waves might propagate along high loops or on 
the surface. It is possible to distinguish this problem usi nR multi-ang;le obse rvations, such as 
the simultaneous observations from SDO and ( STEREO' AKdjisei et al.ll2008l ). Unfortunately, 
the STEREO has not the required high cadence to catch these waves. The results presented 
in Section 3 are based on the assumption that the waves were propagating on the surface. 
Here, we discuss the possibility that the waves were propagating along high loops. Assuming 
the angle of inclination of the high loops relative to the solar surface is 45°, we can obtain 
the propagating average speed of the waves to be ^340 km s~^, and the mean speed of 
the reflected waves ^310 km s~^. These speeds of the waves in the closed transequatorial 
loops are all above the sound speed in the low coronal. By combining the phenomena that 
no signiflcant CME was observed to be associated the waves and that the waves can be 
reflected, we conclude that the waves should be fast magnet osonic waves rather than the so- 
called "pseudo- waves" , which are thought to be the footprints or the low coronal e xtensions 
of the associated CMEs (e.g., IPelannee l[2Q0(^ Ichen et al.lbood : kttriU et al.]l2007l ). 



With the phase speeds of the propagating waves, we can calculate the local magnetic 
field strength of the corona structures where the waves propagated. The sound speed at 
formation temperature of 171 A , i.e., 0.8 MK, is about 131 km s~^, which is calculated from 
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the formula 



147 



[km s ^] (j AschwandenI l2005l ) . with 7 



prup Y IMK 

being the adiabatic index, /x = 1.27 the mean molecular weight in corona (H:He=10:l), 
mp = 1.67 X 10~^^ g the proton mass, and /cb = 1.38 x 10~^^ erg the Boltzmann constant. 
For fast magnetosonic waves propagating along the magnetic field and the amplitude of the 



perturbation perpendicular to the magnetic field, the phase speed is '^ph = 



/47rp 



where 

being the Alfven speed, B the magnetic field strength, and p the ion density. 
Using the phase speed 'Uph = 834 km s~^ of the QFP waves obtained from the k-uj diagram 
of the QF P waves and the num ber density ^ 10^ cm~^ estimated with the 171 A channel 
response (IBoerner et al.ll201ll ). we can estimate the magnetic field strength B of the open 
active region coronal loops which guide the QFP waves, and the result is 5 4.24 Gauss. 
For the waves guided by the closed transequatorial loops, its phase speed '^ph = 240 km s~^ 
is about twice of the sound speed (cg = 131 km s~^). With the same method, we calculate 
the magnetic field strength of the closed transequatorial loops, which is 5 1.24 Gauss, 
about three times smaller than the magnetic field strength of the open active region coronal 
loops. If these waves were propagating in high loops as we discuss above, the magnetic field 
strength of the closed loops should be 1.6 Gauss. 

Previous studies have shown that magnetic reconnection is the basic mechanism respon- 
sible for the fast release of magnetic energy in the corona. However, it is still unclear about 
the origin of the periodicity of solar fiares and the excitation mechanism of the QFP waves. 
Generally, there are two alternative possibiliti es for the g:eneration of the p eriodicity of so- 
lar fiares: spontaneous, and forced periodicity (INakariakov fc Melnikovll2009l ). The former is 
determined by the local plasma properties in the reconnection site, while the latter is excited 
externally. For the first mechanism, magnetic energy is continuously supplied by infiow, and 
it continuously builds up around the vicinity of a fiare epicentre. When the critical point 
for magnetic reconnection is reached, the stored energy is released by a burst, and the same 
process repeats again. A periodic regime of magnetic reconnection can be gen erated in sev- 



eral p hysical situations, such as the coalescence of two magnetic fiux t ubes (ITaiima et al. 
1987 ). the repeated gen eration of plasmoids as well as their coalescen ce ( Kliem et al. 200ol ) 
and other mechanism (lOfman fc Suil l2006l : McLaughlin et al.l l2009l ) . On the other hand 



the periodicity of fiares can a lso be produced by external agents, such as MHD waves (see 
Nakariakov fc Melnikovl l2009l . and references therein). In the case presented here, we con- 
jecture that the periodicity of the fiare might be excited by some internal micro dynamics 
in the reconnection process and it thereby drives the QFP waves as well as the periodicity 
of the accompanying fiare. 

In addition, many authors found evidence that coronal waves are tightly connected wi th 
photospheric global pressure-driven (^-mode) oscillations ( De Moortel et al. 2000 . 2002al lbl: 



Marsh et all I2OO3I : iMarsh fc WalshI bood : Ichen fc PriestI bood : IWang et al.l l2009al lbl).FQr p- 

mode oscillations in the photosphere, they are normally evanescent because their periods are 
well above the cutoff period in the upper photosphere a nd the low chromospher e as well as 
the transition region in nonmagnetic solar atmosphere ( lErdelyi fc FedunI 120071 ) . However, 
for magnetosonic gravity waves, the photospheric p-mode oscillations can leak sufficiently 



into the upper solar atmosphere layers, in which low frequency wave s are channe . 



the highly inclined magnetic fields that can increase the cut off periods (IBel fc Lerov 



This effect has been c onfirmed by the numerical simulation (IDe Pontieul 



observational studi es dsioomfield et a 



Wang et al.lboOQal lbl: IPidkovskv et all 



2006 



Mcintosh fc Jefferies 



200' 



I 



2004a 



ed b y 



bl. l2005h 



Vecchio et al. 



923) 



and 



20071 : 



201 if ). In the present case a few low frequencies are 



revealed by the k-u) diagram of the QFP waves, which do not coincide with the frequencies 
of the pulsation fiare. Here we emphasize that the leakage of p-modes oscillations from the 
photosphere into the low corona should be a noticeable mechanism for driving QFP waves. 



Following IVaughanI ( 120051 ). iGruber et al.l ( I2OIII ) developed a method to test the relia- 
bility of the periodicity of several solar fiares, in which both raw and detrending light curves 
were used to create the power spectra density (PSD) spectrum. They found that the periods 
revealed by the PSD spectrum made from the detrending light curve did not appear on the 
PSD spectrum obtained from the raw light curve, but they did reproduce the instrumental 
periodicity on PSD spectrum obtained from both raw and detrending light curves. They 
concluded that such result should be due to the under-estimating the red-noise component 
in PSD spectrum generated from the detrending light curves, which can introduce spurious 
oscillations by enhancing power of weak signals, we also use this method to test the reliability 
of the periodicity of the fiare analyzed in this paper. We find that the PSD spectrums made 
from the raw light curves do show the similar periods reveale d by the wavelet po wer spec- 
trum, but some of them are lower than the 3a confidence limit ( IGruber et al.ll201ll ). In spite 
of this, we still believe that the periods revealed by the wavelet power spectrums, since they 
all lasted for several distinct cycles during the rising phase of the fiare (see Table [1]) . If the 
periods of the fiare are caused by the erratic, aperiodic red-noise, they should be unstable. 
Finally, we expect to analyze more similar events to confirm the relationship between fiares 
and QFP waves, and theoretical work will be essential to resolve the underlying physical 
mechanism. 
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Fig. 1. — An overview of the wave event, (a) AIA 171 A (b) the filtered image of (a), which is 
obtained by subtracting a smoothed image with a boxcar average over 15 x 15 pixels, (c)-(d) 
AIA 171 A running difference images (Animation 1), and (e)-(f) AIA 193 A base difference 
images (Animation 2). Arrows in frame (a) point to the open and closed transequatorial 
loops which guide the propagating waves, while jointed black arrow in frame (e) indicate the 
multiple arc-shaped wave trains. The white long arrow in frame (c) indicates the direction of 
the associated CME. The white dashed box in frame (e) indicates the selected region where 
Fourier analysis is applied. The coronal dimming regions behind the waves are also indicated 
in frame (f). Sectors A0-A4 in frame (c) and B0-B7 are used to obtain the time-distance 
diagrams along the propagating directions of the propagating waves, and the time-distance 
diagrams are shown in Figure [2] and 7, respectively. The white curve in each frame marks 
the limb of the solar disk (the same in the following figures). The FOV is 590^' x 580^' for 
each frame. 
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Fig. 2. — Time- distance diagrams and analysis of the QFP waves along the open coronal 
loops, (a)-(c) Running difference time-distance diagrams obtained from AIA 171 A running 
difference images along the sectors A1-A3 shown in Figure [I^c), (d)-(f) the same as (a)-(c), 
but they are obtained from the AIA 171 A base difference images. The white arrows indicate 
the QFP wave trains in the time-distance diagram. 
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Fig. 3. — (a) Subregion of Figure [2](b) at times and distances marked by the two red "x" 
signs. The horizontal direction is magnified by 12 times, (b) the intensity profile of each 
pixel (x-axis) as a function of distance (^-axis) at seven consecutive times from 10:56:30 UT 
(red curves). Each profile is normalized into [-6,6] DN and then it is incrementally shifted 
by 12 DN along the x-axis, which equals to the time cadence of the AIA 171 A filter and 
thereby the x-axis also serves as time axis. The blue curves are the fitted curves of a sine 
function, (c) the dashed box region shown in Figure [2](e) , which is also magnified in the 
X-axis direction to fit the panel, (d) the intensity profile of each pixel (y-axis) as a function 
of time (x-axis). Note that the intensity profile of larger distance is stacked on the previous 
one. Each curve is shifted with 5 DN in y-axis direction. 
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Fig. 4. — Fourier analysis of the QFP waves guided by the open loops, (a) Fourier power {k-uj 
diagram) of a three-dimensional data tube of 171 A running difference images during 10:48^ 
11:00 UT in the dashed box region shown in Figure [T](d). The red curve is the normalized 
intensity profile of the straight ridge. The white dotted lines indicate the possible frequencies 
of the QFP waves, which are determined by the peaks showing on the intensity profile of the 
straight ridge, (b)-(d) AIA 171 A running difference sample images in the dashed box region 
shown in Figure [T](d), they are Fourier filtered with a narrow Gaussian function centered 
at the peak frequencies of 11.2 mHz, 15.6 mHz, and 23.7 mHz respectively. The FOV for 
frames (b)-(d) is 200'' x 200''. 




Fig. 5. — Sample light curves of the accompanying flare, (a) time proflle of GOES 1-8 A 
SXR. (b) RHESSI HXR count rates in the energy band (4 seconds integration) of 12-25 
keV (green), and AIA 171 A (purple) and 1600 A (red) light curves measured from the flare 
region. The time range is from 10:48-11:00 UT. (c) the same with (b), but for detrended 
fluxes obtained by subtracting the smoothed fluxes using a 120 s boxcar. 
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Fig. 6. — Wavelet analysis of the detrended light curves of the flare in UV, EUV, and hard X- 
ray. For each wavelet power spectrum, the corresponding normalized global power is plotted 
on the right, in which the dotted line indicates the 95% signiflcance level. The red horizontal 
dashed lines indicate the possible periods of the flare, which are determined by the peaks 
(above 95% signiflcance level) of the global power. The blue contours in each wavelet power 
spectrum outline the region where the signiflcance level is above 95%. In this flgure, redder 
color correspond to higher wavelet power. 
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Fig. 7. — Time-distance diagrams for analyzing the waves guided by the closed transequato- 
rial loops, (a)-(d) AIA 193 A running difference time-distance diagrams obtained from the 
sectors B2-B5 shown in Figure [I](e) . The white horizontal and the oblique arrows point to 
the reflected wave trains, while the jointed white arrows in frame (b) indicate the multiple 
propagating wave trains. 
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Table 1: All the periods, frequencies, duration, cycles, averaged periods and frequencies of 
the flare 



Wavelength 


Pf (s) 


7~i / TT \ 

Ff (mHz) 


Duration 


1 

Cycles 


-Pfa (s) 


7~i / TT \ 

Ffa (mHz) 


A T A -1 >~7 r\ r\ A 

AIA 1700 A 


17 ± 4 


59 


264 


15.9 


17 ± 4 


59 


A T A ^ ^ ^ A 

AIA 171 A 


25 ± 2 


40 


170 


6.8 






AIA 211 A 


25 ± 3 


40 


166 


6.6 


25 ± 2.5 


40 


AIA 171 A 


31 ± 3 


32 


288 


7.4 






AIA 193 A 


32 ± 4 


31 


219 


6.8 






A X A £^ c\ ^ A 

AIA 335 A 


32 ± 2 


31 


160 


5.0 






A X A £^ r\ A A 

AIA 304 A 


33 ± 4 


30 


246 


7.5 






A X A -1 r\ r\ A 

AIA 1700 A 


35 ± 2 


29 


303 


8.7 


32.6 ± 3 


31 


AIA 193 A 


40 ± 5 


25 


269 


6.7 






A T A r\ A A 

AIA 304 A 


40 ± 2 


25 


246 


6.2 






A X A -1 A 

AIA 1600 A 


41 ± 5 


24 


269 


6.6 






A T A t l~r r\ r\ A 

AIA 1700 A 


42 ± 3 


24 


282 


6.7 


41 ± 3.8 


24.4 


A X A ^ l~r ^ A 

AIA 171 A 


53 ± 4 


19 


296 


5.6 






AIA 304 A 


53 ±5 


19 


291 


5.5 






AIA 335 A 


53 ±5 


19 


271 


5.1 






AIA 193 A 


54 ±9 


19 


314 


5.8 






AIA 211 A 


54 ±5 


19 


305 


5.6 






AIA 94 A 


56 ±5 


18 


303 


5.4 






AIA 131 A 


56 ±5 


18 


299 


5.3 






AIA 1700 A 


56 ±8 


18 


303 


5.4 






AIA 1600 A 


57 ±7 


18 


306 


5.4 






RHESSI 12-25 keV 


62 ±5 


16 


307 


5.0 






AIA 1700 A 


80± 15 


13 


330 


4.1 


55.4 ±5.8 


18 


RHESSI 12-25 keV 


86± 12 


12 


325 


3.8 


83 ± 13.5 


12 


RHESSI 12-25 keV 


143 ± 19 


7 


327 


2.3 


143 ± 19 


7 



Note. — In this table, Pf and Ff are the periods and frequencies of the flare in various wavelength bands, 
while Pfa and Pfa are the averaged periods and frequencies respectively. 



